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ABSTRACT 


Formulas  are  derived  for  volume  scattering  coeffi¬ 
cients  applicable  to  aerosols.  Different  types  of  size 
distribution  for  the  aerosol,  such  as  normal  distribu¬ 
tion,  log-normal  distribution,  and  exponential  distri¬ 
bution,  are  assumed. 

Numerical  results  obtained  by  the  author  and  those  pub 
lished  in  the  literature  are  given  and  discussed. 

The  volume  scattering  coefficients  become  smooth 
functions  if  the  size  distribution  is  very  wide.  Hence, 
simpler  mathematical  formulations  can  be  used  than 
those  obtained  from  the  exact  Mie  theory.  However, 
only  the  number  of  aerosol  particles  can  be  obtained 
from  an  evaluation  of  experimental  data,  but  details 
of  the  size  distribution  cannot  be  extracted  from  ex¬ 
perimental  data. 
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i.  INTRODUCTION 


Thr  scattering  functions  for  single  particles  with  real  index  of  refraction  n  are 
known  for  many  values  ofn  ,  and  a  very  large  range  of  size  parameters 
a{a  s  2rrt/\  ).  Total  scattering  coefficient,  angular  scattering  coefficients,  and 
degree  of  polarization  are  known.  Because  of  the  electronic  computer  technique, 
such  computations  pose  no  problem  and  scattering  coefficients  can  be  generated 
for  any  desired  value  of  n  and  size  range.  For  aerosol  studies,  where  integra¬ 
tions  over  a  particle  size  range  and  wavelength  range  are  needed,  the  minor 
fluctuations,  the  so-called  ripples,  are  unimportant.  Smoothed  values  without 
ripples  are  often  completely  sufficient.  Hence,  this  problem  is  well  understood 
(van  deHulst,  1957). 

The  scattering,  absorption,  and  polarization  coefficients  for  single  particles 
with  complex  index  of  refraction  “rr  =  n  -i*  are  also  known  for  many  values  of 
n'anda  .  However,  less  extensive  data  have  been  generated  than  for  real  n. 

The  results  are  well  understood;  absorption  reduces  the  fluctuations,  the  ripples 
vanish  first,  and  the  major  fluctuations  become  weaker  the  larger  the  absorp¬ 
tion,  i.  e. ,  k  .  There  exists  a  steep  increase  of  the  scattering  coefficients  in 
the  Rayleigh  range  for  «  <  1  and  a  very  simple  functional  behavior  for  a  >  1. 
This  applies  to  the  total  as  well  as  to  the  angular  scattering  coefficients.  Again 
the  problem  is  understood;  data  exist  to  derive  the  general  behavior  and  machine 
computations  pose  no  difficulties  (Deirmendjian,  1963). 

The  most  important  problem  of  today  is  the  scattering  from  a  small  volume 
containing  an  aerosol  of  uniform  refractive  index  o'  or  n  .  The  particles  inside 
the  volume  are  not  of  uniform  size,  but  possess  a  size  distribution.  The  size 
distribution  can  be  described  by  reasonably  simple  functions.  As  such,  one  may 
assume  an  exponential  law  (as  in  Junge's  distribution  function),  or  a  Gaussian 
distribution  function,  or  any  other  which  has  been  proven  to  exist  (Deirmendjian, 
1963;  Dettmar  et  al,  1963),  A  literature  survey  was  published  by  Penndorf 
(1963). 

In  this  report  we  discuss  some  theoretical  derivations  for  representing  the 
scattering  coefficients  for  assumed  size  distributions  and  show  some  of  the 
results  of  such  computations. 


2.  total  and  angular  volume -scattering  coefficients 

The  total  volume  scattering  coefficient  f}j  for  an  aerosol  with  uniform  size 
distribution  is 

0T  =  nt2  K  l  a,  n  I  N  U) 

where 

r  =  radius  of  aerosol  particle 
a  -  size  parameter;  a  =  ( 2nt)/\  =  kr 
N  =  number  of  particles  of  radius  r 
n  =  refractive  index  of  aerosol  particle 
K  =  total  Mie  scattering  coefficient. 

The  angular  volume  scattering  coefficient  is  correspondingly 

00  -  n  i0  la,  n,  0\  N  (2) 


where  is  the  angular  Mie  scattering  coefficient  for  a  single  particle  in  the 
direction  0;  0  =  0  degree  for  forward  scattering. 

The  relation  between  ig  and  K  is  (see  Penndorf,  I960) 

4  n 

i  0  d  w  ■  K  .  (3) 


Hence 


(4) 


In  the  following,  the  formulas  are  written  for  K  but  by  replacing  K  with  i0  the 
results  for  0T  follow  automatically  for  ft#  \  we  used  the  symbol  0  for  0T  in  the 
following  for  simplification  only. 


The  computations  of  the  volume -scattering  coefficient  for  an  aerosol  containing 
particles  of  non-uniform  sires;  i.  e.  ,  and  aerosol  with  a  specific  size  distribution 
follows  from 

d/S  7  dN; 

- —  =  >/r?  K  |aj,  n  |  —■ (5) 
da  da 

1  1 

where  rj  =  radius  at  r  ±  ( —  )  Ar  or  a  ±  (—  )  Aa ,  and  dN;/da  =  number  of 
particles  between  a-t  ±  A  a. 

This  can  be  written  as 

d/3  ,  1  dN  na2  dN 

da  i  .  kj*  da  k*  da 

It  is  understood  that  all  these  consideratio  is  refer  to  a  specific  wavelength  A, 
and  the  integration  over  the  selected  wavelength  range  should  be  performed  if 
a  wide  band  is  chosen.  We  will  come  back  to  it  later.  It  is  also  understood 
that  n  remains  unchanged  for  the  full  wavelength  range  selected. 

We  write  the  size  distribution  function  p  I  I  as 


=  p  lr  I 


or  with  dr  =-j-  da 
dN  1 

*r~  "T  p  1*1  * 

da  k 


The  total  number  of  particles  follows  from  equation  (8)  as 


N  .  /dN  .1  J 


p  |a|  da  . 


For  p  |a|  we  shall  later  on  introduce  specific  population  functions. 


Using  equation  (8)  we  can  now  write 
dp  "a2  .  . 

_  Kp  laj  . 


a  relationship  which  holds  for  a  fixed  interval  da  *nd  a  fixed  wavelength  * . 


The  integration  over  all  "a yields 
/3  = -  f  a2  K  la|  p  |a|  da  . 


(ID 


In  this  form  /3  is  valid  for  a  specific  X  only,  because  k  contains  X. 

Next  we  define  a  new  total  Mie  coefficient  Kp  for  a  size  distribution,  so  that 
equation  (1)  can  be  used 


it  a 


P  =  —  Kp  K(N  • 

k2 


This  equation  (12)  is  identical  to  equation  (11)  and  it  follows 


(12) 


h  K'  •  - 


kN  . 

0  0 


/ 


a2  K  |a|  p  |a|  da 


(13) 


In  this  form  aQ  denotes  the  a  value  at  which  the  distribution  function  reaches  a 
maximum;  i.  e.  ,  dN/da  reaches  a  maximum. 

Introducing  n  from  equation  (9)  into  equation  (13)  leads  to 


Kp  K* 

a* 

o 


/  a2  K*al  p  la!  da 

4 


/ 


(14) 


p  |a|  da 


In  the  case  of  angular  scattering  consideration  we  define  an  angular  Mie 
scattering  coefficient  as 


/ 


a2  I#  tat  p  )a|  da 


*0p 


/ 


(15) 


p  (aide 
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3.  AVERAGE  SIZE 


3.  1  AVERAGE  SIZE  a  FOR  p  la} 

The  average  radius  r  of  particles  having  a  size  distribution  p  Irl  is  the  arithme 
tic  mean  value  of  all  the  radii,  if  N  is  the  total  number  in  the  sample. 


Thus 


r 


(fj  +t2  +  ... 


(16) 


Replacing  r  by  o^k  and  the  sum  by  the  integral  leads  to 
ap  la} da  . 


a  a 


-  f 

kN  / 


(17) 


Since 


/ 


p  |a|  da  , 


ap  latda 


a  • 


/ 


p la} da 


3.  2  AVERAGE  SIZE  a  FOR  NORMAL  DISTRIBUTION 
The  normal  distribution  is 


» ^  •*  (-  if) 


with 


OB) 


09) 


v  •  (•  -  aj/l  »d  a  •  •»  ♦  % 


Equation  (3  8)  can  be  solved 


and 


ap  1  a  I  d  c 


p  |  a\  da  =  —  . 

a 


+ 


Hence 


a 


%  +  w 


+  0.798  ff 


(20) 


Thus,  for  a  normal  distribution,  the  average  size  is  determined  by  the  and 
the  standard  deviation  a,  a  is  always  larger  than  %  . 


( 


4.  Kp  AND  iflp  FOR  NORMAL  DISTRIBUTION 


We  select  a  population  function  for  which  the  aerosol  particles  have  a  normal 
distribution  with  respect  to  a.  The  total  number  of  particles  in  the  sample  is  N 


P  I  a  I 


N 

OypZ 


e*p  (-u2/2«)  m  (N/a)  $  Ju  | 


(21) 


whe^e  ^iu  1  denotes  the  normal  distribution. 

The  maximum  of  the  distribution  is  given  by 

p  l  -  (N/«r/2»)  »  (0.399  N /a)  .  ;  (22) 

Equation  (14),  using  the  result  of  the  last  section,  namely 


becomes 


/ 


Kp  1  a  j  *  - —  j  a2  K  1  a  1  p  1  a  !  da 
-0 


cr 


kpK' 


7  /  °2ki°i*{t(!71)  }d° 
0  0 


(23) 


The  integral  has  to  be  solved  numerically  because  K  is  not  a  simple  function  of 
a;  ^(u)  is  tabulated  in  standard  text  books  on  statistics. 

In  the  case  of  angular  scattering  consideration,  Kp  is  replaced  in  equation  (23) 
by  i  ^  and  K  by  .  The  general  formalism  remains  unaltered. 

In  many  aerosols  the  so  called  logarithmic  distribution  (log-normal)  has  been 
found.  In  this  case  the  q’lantity  u  in  equation  (21)  is  defined  as 


or 


u  «  [ 


a  m  ^  exp  (ou) 


(24) 


We  can  shelve 


y*  eo 

J  p  l«  Id  a 


which  is 


f 


P  1 « I 


I  ,2) 

0 


[ln(a/aoy]  j  da 


(25) 


/M  2 

em  e“Am 

_ 

is  given  by  Grobner  and  Hofreiter  (1950,  p.  57) 

^  p  le|  ■  N  aQ  |cxp  ^  [l  -  <f>  o/\J~2)] 


dm  for  which  a  solution 


(26) 


/*  2 

e-*  ,  the  error  intregai,  which  is  tabulated  in  standard 

textbooks. 

The  solution  for  angular  scattering  is  now  obtained  from  equations  (15)  and 
(26)  as 


i  q  |  g  |  »  '  -  -  -  — 

°  V2  (**P  (» V 2  )  1  [  1  -  i>  (—  a/\f2 )  ] 


(27) 


Because  of  the  Mie  function  in  ,  a  numerical  evaluation  of  the  integral  has 
to  be  performed. 


da 
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5.  EXPONENTIAL  POPULATION  FUNCTION 


Junge  (1963)  has  shown  that  a  good  agreement  for  the  natural  aerosol  in  the 
range  0.  1  <  r  <  20p  is  given  by 

_fL_  .  N  ry  (28! 

d  dog  r)  ° 


A  good  agreement  is  normally  found  for  y  = 

\y\  >  IU  • 


-3.  In  all  cases  y  is  negative  and 


Such  a  function  leads  to 


-  -  (log  e)  r^  *  . 

dr  0  * 


Replacing  r  by  o 


dN  / 1  \  y  v  ,  „  , 

P  -  —  -  Nq  (log  e)  f- j  a>  -1  =  CN0  aY~l  . 


Such  a  exponential  size  distribution  is  only  valid  within  specified  limits  of  r 
(and  therefore  a).  Hence, 


a^  K  laj  p  |a{  da 


1  U1 

KP  '“o'  -  - - - 

”o  t  2 


p  1  a |  da 


The  denominator  can  be  written  as 


p  la|  da  =  CN 


(a/  -  a,/)  • 


Since  the  number  of  particles  at  the  upper  end  (large  r )  is  very  small  and 
determined  by  the  limitations  of  the  instrument,  we  can  set  the  integral  equal 
to 


-9- 


Thus  equation  (30)  becomes 


a  ^  K  |al  da  . 


(35) 


and  for  y  =  -4 


Kplal‘  =  4«12  j 


Similar  formulas  can  be  written  for 


;  for  example,  for  y  =  -3  one  obtains 


(36) 


It  follows  from  the  foregoing  theoretical  considerations  that  we  can  compute 
Mie  scattering  functions  for  aerosols,  provided  the  general  size  distribution 
law  can  be  assumed  from  other  observations  or  general  physical  considerations 
If,  however,  no  knowledge  on  the  general  size  distribution  law  exists,  the 
experimental  observations  cannot  be  interpreted. 
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6.  RESULTS  OF  COMPUTATIONS 


Using  the  theoretical  derivations,  we  show  in  several  examples  what  values 
are  obtained  for  Kp  and  igp . 

For  a  normal  distribution,  equation  (23),  we  choose  =  5,  <*  =  1.  0,  and 
da  =  0.  2,  and  use  the  K  values  for  a  =  1.  33  obtained  by  Penndorf  (1956). 

In  this  case  we  obtain  Kp  !  5.  0  I  =  3.  69.  The  value  for  KI5.  0  l  =  3.  59.  In  table  I 
we  give  some  examples  for  n  =  J.  33  and  aQ  -  6  0(0.  2)  7.  0. 

TABLE  I 

TOTAL  SCATTERING  COEFFICIENTS  F_  FOR 
SOME  SIZE  DISTRIBUTIONS  P 
WITH  a  =  1.0 


ao  6-  0 

6.  2 

6.  4 

6.6 

6.  8 

7.  0 

Kp  3.  84 

3.  82 

3.  83 

3.  75 

3.70 

3.63 

K  3.89 

3.89 

3.  96 

3.  96 

3.84 

3.74 

As  expected,  smooth  variations  in  Kp  result  for  size  distributions. 

Kerker  (1964)  conducted  experiments  and  theoretical  computations  for  refrac¬ 
tive  indices  (n  =  1.  43,  1.  51,  and  2.  07),  a  size  range  «0  =  1.  9  (0,  1)  15,  and 
width  parameter  a o  =  0.  000  (0.  05)  0.  155. 

He  used  a  logarithmic  normal  distribution  and  aQ  and  crQ  are  the  parameter  for 
such  a  size  distribution. 

Experimental  data  are  obtained  and  compared  with  the  theoretical  computations. 
He  found  that  the  experimental  angular  scattering  data  for  the  two  polarized 
components  in  the  range  6  =  30  to  130  degrees  can  be  fitted  to  the  theoretical 
data;  he  obtained  the  modal  size  and  the  width  parameter  of  the  distribution 
function.  The  uniqueness  of  the  solution  has  not  posed  a  problem.  The 
experimental  data  always  fitted  one  set  of  theoretical  data  with  the  least  devia¬ 
tions. 

His  results  show  that  for  very  narrow  distribution  (aQ  <  0.  155)  theoretical 
computations  can  be  applied  to  experimental  data  and  the  important  parameters; 
namely,  refractive  index  n  and  size -distribution  parameters  will  be  obtained.  - 

For  atmospheric  aerosols  the  situation  is  not  so  good,  because  we  have  a 
mixture  of  refractive  indices  and  very  broad  size  distribution. 
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For  an  exponential  distribution  (Junge's  aerosol  distribution)  we  have  computed 
data  based  on  our  theoretical  results  and  find  the  following  values  for  Kp  {see 
table  II)  for  the  interval  from  =  0.  1  to  =  50  using  -  -3, 

TABLE  II 


FOR  EXPONENTIAL  DISTRIBUTION 
WITH  Kp=-3 


n 

1.  33 

1.  40 

1.44 

1.  486 

1.  5 

2.  0 

kp 

1.  26 

i.  48 

1.61 

1.  75 

1.  79 

3.  13 

The  value  for  n  =  2.  0  is  not  so  accurate  as  those  for  the  interval  n  =  1.  33  to 
1.5. 

The  effect  of  different  values  of  y ,  i.  e.  ,  different  steepness  of  the  size  dis¬ 
tribution  has  been  computed  for  the  internal  a  =  0.  1  to  a  =  30  and  n  =  1.  33. 
We  find  for  y  -  -2,  Kp  =  1.  93;  y  =  -3,  Kp  =  1.  17,  a.id  for  y  =  -4,  Kp  =  0.  38. 

Such  a  result  is  expected  because  the  steeper  the  function  the  more  aerosol 
particles  are  contained  in  the  very  small  a  range. 
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